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We develop a non-perturbative method that yields analytical expressions for the deflection angle 
of light in a general static and spherically symmetric metric. The method works by introducing 
in the problem an artificial parameter, called 5, and by performing an expansion in this parameter 
to a given order. The results obtained are analytical and non-perturbative because they do not 
correspond to a polynomial expression in the physical parameters. Already to first order in 5 the 
analytical formulas obtained using our method provide at the same time accurate approximations 
both at large distances (weak defiection limit) and at distances close to the photon sphere (strong 
deflection limit). We have applied our technique to different metrics and verified that the error is at 
most 0.5% for all regimes. We have also proposed an alternative approach which provides simpler 
formulas, although with larger errors. 
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I. INTRODUCTION 

The theory of General Relativity (GR) predicts that massive bodies deform the spaee-time around them; even 
massless particles, like photons, will therefore feel the gravitational force and their trajectories will necessarily depart 
from a straight line. Such effect is particularly strong in the proximity of black holes, which are very massive and 
compact objects: in fact at a certain distance from the black hole (the photon sphere) the deflection angle of the 
photon becomes infinite. This regime is known as strong deflection limit (SDL), whereas the regime corresponding to 
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deflection at large distances is known as weak deflection limit (WDL). 

The WDL has been studied in a series of articles for different metrics: for example the lensine; from a Kerr metric 
has been considered in whereas the lensing from a Reissner-Nordstrom metric has been considered in 

In recent times there has also been great interest in studying the effects of strong deflection limit, i.e. lensing due 
to light passing very close to a compact and massive body. For example SDL in a Schwarzschild black hole has been 
considered by Frittelli. Kling and Newman Q and by Virbhadra and Ellis Q; Virbhadra and collaborators have also 
treated the SDL by naked singularities and in the presence of a scalar field Eiroa, Romero and Torres [lol | 
have described Reissner-Nordstrom black hole lensing, while Bhadra has considered the gravitational lensing du e 
to the GMGHS charged black hole 



11 



^ ^, Bozza has studied the gravitational lensing by a spinning black hole [l^ : 

Whisker 1^ and Eiroa [3| have considered SDL by a brancworld black hole; still Eiroa [l^ has recently considered 
the gravitational lensing by an Einstcin-Born-Infeld black hole; Sarkar and Bhadra have studied the SDL in the 
Brans-Dicke theory [l^ ; Konoplya has studied the corrections to the deflection angle and time delay for a black hole 
immersed in a uniform magnetic field l3|; Gyulchev and Yazadjiev have studied the SDL for a Kcrr-Sen dilaton axion 



black hole [18[; finally Perlick 



19| has obtained an exact gravitational lens equation in a spherically symmetric and 



static spacetime and used it to study lensing by a Barriola-Vilenkin monopole and by an Ellis wormhole. Notice that 
Bozza and Scrcno [20| have also investigated the SDL of gravitational lensing by a Schwarzschild black hole embedded 



in an external gravitational field. 

We can distinguish between two different approaches that have been developed to obtain analytical expressions 
for the deflection angle in the strong regime: one which looks for improvements of the weak lensing expressions, 
whose range of validity is therefore extendend to distances closer to the photon sphere, without however taking into 
account the divergence of the deflection angle on the photon sphere, and one which treats exactly the singularity of 
the photon sphere and whose precision rapidly drops at larger distances ^. In the flrst category falls the work of 
Mutka and Mahonen 221 and of Belorobodov 23| who worked out improved formulas for the deflection angle in a 
Schwarzschild metric, and the more systematic approach of Keeton and Fetters ^ who have developed a formalism 
for computing corrections to lensing observables in a static and spherically symmetric metric beyond the WDL. In 
the second category falls the work of Bozza, who has introduced an analytical method based on a careful description 
of the logarithmic divergence of the deflection angle which allows one to discriminate among different types of black 
holes 24 1 . Recently, Iyer and Fetters 2^ have also developed an analytic perturbation framework for calculating the 



bending angle of light rays traversing the fleld of a Schwarzchild black hole, obtaining accurate expressions even in 
proximity of the photon sphere. 

In a different category falls a method developed by Amorc and collaborators 
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27[. This method enables one to 



convert the integral for the deffection angle in a static and spherically symmetric metric into a series in an artiflcial 
parameter i5. Such series has an exponential convergence rate and its terms can be calculated analytically. The method 
is non-perturbative in the sense that it provides non-polynomial expressions in terms of the chosen physical parameter 
and yields sufficiently accurate results even at first order. In our previous works we have tested our formalism on a 
variety of different metrics, always obtaining very encouraging results. 



^ For a detailed discussions on the photon surface the reader can refer to [21 
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The purpose of this paper is to improve our method in order to provide an accurate treatment close to the photon 
sphere, even at first order, but without the customary deterioration of the resuhs at larger distances. 

The paper is organized as follows: in Section [IT] we illustrate the application of our method by means of the 
Schwarzschild metric and later we show how to treat a general case. In Section Hill we extend our analysis to a general 
static and spherically symmetric metric. In Section llVI we compare our approximations with exact results and other 
approaches available in the literature. Finally, in the last section we briefly summarize and discuss our results and 
consider further developments. 



II. FORMALISM 



Let us first review the method of Amore and collaborators 
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27|. We are interested in the general static and 



spherically symmetric metric which corresponds to the line element (in the following we set the velocity of light c = 1) 

= B{r)dt'^ - A{r)dr'^ - D{ry [dO"^ + sin^ 6* dcf,'^) (1) 

and which contains the Schwarzschild metric as a special case. We also assume that the flat spacetime is recovered 
at infinity, i.e. that limr^oo /(?') = 1, where /(r) ~ {A{r), B{t), D{r)) ^. 

The angle of deflection of light propagating in this metric can be expressed by means of the integral ^ 



A(j> 



^A{r)/D{r) 



\ 



D{r) Bjro) 
roj D{ro) B{r) 



- 1 



dr 
r 



(2) 



where tq is the distance of closest approach of the light to the center of the gravitational attraction. 
By introducing the variable z = r^/r one can rewrite the equation for the deflection angle as 

dz 



where 



V{z) 



/o ^V{l)-V{z) 
^Dir^lz) D^iro/z)B{ro) 



71", 



B{ro) 



(3) 



(4) 



A{ro/z) Airo/z)B{ro/z)D{ro) D{ra) 
is a sort of "potential" built out of the metric. Notice that the integral in Eq. ([3]) can be solved analyticallyonly in a 
limited number of cases, such as for the Schwarzschild metric jsil l and for the Reissner-Nordtrom metric where 
it can be expressed in terms of elliptic integrals. No analytical formula can can be obtained in the case of a general 
static and spherically symmetric metric. 

Since an explicit calculation of the integral is not possible, we interpolate the actual potential V{z) with a simpler 
potential Vo{z), which should be chosen in such a way that the integral ([3]) can be performed explicitly when V{z) is 
replaced with Vq{z). Then we write 



Vs{z) = Vo{z)+d{V{z)~Voiz)), 



(5) 



^ This is a sufRcient but not necessary condition since it warrants that /(r) is analytic around r = oo; for example in [29| we have applied 
the method to the Weyl metric which is not asymptotically flat. 
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where 5 is a dummy parameter. In general Vo{z) may depend upon one or more arbitrary parameters; for the time 
being we simply choose Va{z) — Xz^. We can rewrite the expression for the deflection angle as 

AcPs = 2 /' , , ^ - TT, (6) 

Jo ^Vo{l)~Vo{z) v/l + '5A(z) 

where 

Notice that Eq. ^ reduces to Eq. ^ ior 5 — 1 and therefore is not an approximation. The expansion of EqJ6]) in 



powers of S converges at (5 = 1 provided that |A(z)| < 1 for < z < 1. As discussed in earlier papers j26l. l28l. l29l| this 
condition requires that A be greater than a critical value Ac- In that case one obtains a parameter-dependent series 
that converges towards the exact result which is however independent of A. The artificial dependence on A observed 
in the partial sums A0^^\ N ~ 1,2,..., is minimized by means of the Principle of Minimal Sensitivity (PMS) [s^, 
which corresponds to imposing the condition 

^A^(^) = 0. (8) 



26| 



A proof of convergence of the series and an estimate of its rate of convergence are given elsewhere 

One might be tempted at this point to question our definition of the method as being non-perturbative: after all, 
the method works by performing a perturbative expansion in S. However, one should understand that the solution of 
Eq. ([8]) is in general a function of the parameters in the problem and when substituted back in the series it provides 
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non-polynomial expressions in the physical parameters 
it is set to one at the end of the calculation. 

We are now ready to generalize this method. The first step is to write the integral as 

r dz „ dz 

A0 = 2 / , + 2 



27l |. whereas the dependence upon 5 desappears because 



^V{l)~V{z) Ja ^V{l)~V{z) 
= AcPa + A(Pb-n, (9) 

where cr g (0,1) is an arbitrary point in the region of integration. The two integrals in this equation will now be 
approximated following two different strategies. Clearly, the particular case a = I corresponds to the method just 
outlined 



26 



27| 



For clarity, we confine ourselves, for the moment being, to the case of the Schwarzschild metric and later we 
generalize our results to arbitrary metrics. The Schwarzschild metric is given by 



and the potential V{z) reads 



B{r)^A-'{r)=[l-^Y D{r) ^ 1, (10) 



V{z) = z'-—z\ (11) 
3^ 



where /i = ro/3GM > 1. Fig. [T] shows the potential V{z) with /i = 3/2. The dashed and dotted lines correspond 
to the quadratic Taylor polynomials around z = and z ~ fi. Viewed from the perspective of a classical mechanics 
problem, the points z ~ and z — j.i correspond to a stable and an unstable point of equilibrium, respectively; the 
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FIG. 1: The potential V{z) in Eq. with /i = 3/2 (solid line). The dashed and dotted lines correspond to quadratic Taylor 
expansions around z — and z — ^, respectively, (color online) 



nonlinear pendulum, for example, is a simple physical system that displays this behavior. Since the integral that we 
want to calculate is restricted to z < 1, the unstable point of equilibrium will not be reached unless fi = 1. In such a 
case the integral will diverge and it will correspond to the photon sphere'^. 



One expects ct « 1 to be a reasonable choice when /i ^ 1, since it has given accurate results earlier [26|, |27[. On 
the other hand, as /.t 1+ the second integral will become increasingly more important and one expects the optimal 
value of <T to move to the center of the integration region. We will see that this is the case later on when we discuss 
a systematic way of partitioning the integral. 

Our strategy is simple and consists of treating the integral in each region differently. In the first region we follow 



essentially the earlier procedure 26|, |27[ with the interpolating potential Va{z) = Az^, A > 0. In the second region, on 
the other hand, we will interpolate the potential with the inverted parabola Vf,(z) = V{^) — p{z — /i)^, where p > 
is another arbitrary parameter. Notice therefore that we are working with three arbitrary parameters, A and p which 
enter in the definition of the interpolating potentials, and a which defines the border between the two regions. 
In the first region we write 

dz 1 



where 



Aa(z) 



V{l)~V{z) _ 3(A-l)Ai(z + l) + 2(z2 + z+l) 



K(l)-K(^) 
After expanding to first order in 6 we obtain 



dz 



1 - 



3Xp{z + l) 



(12) 



(13) 



(14) 



^ Notice that fJ, = 1 yields the well-known result ro = 3GM. 
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Straightforward integration yields 



^) _ -2a (VI - CT^ - 2) + 3(3A - l)/i(cr + 1) arcsin(cr) - 4Vl - cr^ + 4 



3A3/2^(a + l) 



The PMS (see Eq. ([8])) gives us the optimal value of A 



^PMS 



2a (VI - cr^ - 2) + 4 (Vl - 0-2 - 1) + 3Ai(cr + 1) arcsin(cr) 



and 



3/i(CT + 1) arcsin(cr) 
2 arcsin((T) 



ia + 2)Vl-(T^ -2{a + l] 



\ 3/i(a + l) 



(15) 



(16) 



(17) 



arcsm cr 



In the second region we have 



A<Pb = 2 



where 



^ V(l)- V(z) _ _ -6p/i^ + 3(p + l)(z + 1)m - 2 (z2 + z + 1) 



Vbil)~Vbiz) 3np{2fi-z-l) 



After expanding to first order and integrating we obtain 



2 (e +^i-^ia + a-l)+ 3A*(3p - 1)^ In [ _^^^J\^ 



where we have defined 



C = Vcr2 _ 2//C7 + 2/i- 1. 



(18) 



(19) 



(20) 



(21) 



It is worth noticing that this simple first-order approximation exhibits the correct logarithmic singularity at ^ = 1. 
The PMS gives us again the optimal value of p, 

2{-3^i + a + 2)^ 



Ppms — 1 



3fi{2fi — (7 — 1) In 



1-M 

-/i + (T + C 



(22) 



and 



A0« 



By adding the two expressions we obtain 



2V31n( 

V -/i + g + 4 



\ 



2(-3/i + g + 2)g 



(23) 



^(2^ - (T - l)ln 



2x/3 VmM~T) arcsin^/^ (cr) 



y^2cr (Vl - cr2 - 2) + 4 (Vl - (T^ - l) + 3/i(cr + 1) arcsin(cr) 



+ - 



2V31nf^-i^ 



71", 



(24) 
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which still depends on the arbitrary parameter a. 

Figs. [2] and [3] show the approximate deflection angle given by Eq. ([24|) and its percent error S = 100 x 
(A</)p]^g — A(^exact)/A(^oxact for two valucs of /i. We appreciate that the arbitrary parameter a can also be de- 
termined by the PMS. In Fig.[2l ^ = 10 and the optimal value of a is close to 1, as expected. On the other hand, if we 
take /i = 1.001, i.e. close to the photon sphere, the optimal value of a drops to about 1/2. It is quite remarkable that 
in both cases the error made by choosing the optimal value for g is smaller than 1%. Another important observation 
is that the maximum of A(j) is quite Rat and, consequently, a shghtly imprecise estimation of apMS not affect the 
precision of the approximation drastically. 

For this reason we do not pretend to obtain a directly by solving the PMS condition, dAcp'-^^ /da = (which is 
equivalent to finding the maximum of the curve in the left panels of Fig. [2] and [3]), since that would certainly be 
a difficult task and lead to quite involved expressions, but we rather use a simple analytical approximation, which 
correctly describes the limits fi oo and /i ^ 1^. As noticed above, since the maximum is quite flat for /i ^ 1+ one 
expects only a modest loss in precision, while providing much simpler expressions. 

In Fig. m we have plotted the exact value of cr obtained solving numerically the PMS condition and the reasonable 
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FIG. 4: (TpMS obtained numerically solving the PMS condition (broken line) as a function of l//.t and the linear approximation 
a = 1 — 1/2^ (solid line), (color online) 
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FIG. 5; Percent error of Eq. ([25} as a function of fi for the deflection angle in the case of the Schwarzschild metric. Left panel; 
range 1 < /i < 2; Right panel: range 1 < ^ < 1000. 



analytical approximation crpMS ~ 1 — 1/2/i. Thus we obtain 



arcsin 1 



3/2 



6 arcsin (l ^ ^) - + ^6 In (^) - 6^ + ^ + 3 



TT, 



(25) 



which provides an accuracy better than 1% for all values of fi, even arbitrarily close to the photon sphere, as shown 
in Fig.O Later on we will derive an even simpler analytical formula for the deflection angle from Eq. (|25p . 



III. GENERAL METRIC 



We will now attack the problem of obtaining a first order formula for a general static and spherically symmetric 
metric, in analogy with what has been done in Refs. 26l. l27j. For a given metric, once the functions A{r), B{r) and 
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D{r) are given (sec Eq. ((T|)), one obtains a potential V{z), as previously explained. We assume that the potential 
admits two different expansions, one at 2; 0, and one around z = fi, which is a local maximum: 



Viz) = VnZ'' 



(26) 



Notice that the first series runs from n = 1, since z ~ is not necessarily a local minimum of V{z), although this 
was the case for the Schwarzschild metric. Clearly, in particular cases, such as the one previously examined, the series 
coefficients vanish after a certain value of n, thus yielding polynomial potentials. 

Following the discussion of the preceding section we split the integral as in Eq. © and proceed to the calculation 
of each part. 

1. Region I (0 < z < a) 



In this region we expand V{z) around z = and use the interpolating potential Va{z) ~ Xz^. After expanding to 
first order we obtain Eq. p4|) with 

00 n— 1 ^ 



ri=l fe=0 



The deflection angle can now be written as 



oc n— 1 



A(^i^' = — arcsin (a) - —j^ ^ w„ ^ 4(cr) 



ri=l fc=0 



where we have defined the integrals 



arcsin^ sin*^' I 



f + sin ( 



-d0 



which can be calculated exactly; for example: 



l2{<y 



= I- 



arcsm(a-) H j= — 1 

Vl + 



arcsin((T) — \/ 1 — a'^ + 2 

Vl + cr 



^ arcsin(a) - 2 + ^ {-a^ + a + 4) J - 1 



I l-a 



- arcsin(fT) H h -W (-2a-^ + - 7a - 

2 ' 3 6 V CT + 1 ^ 
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The PMS yields 



so that 



^PMS — 



arcsin (a) 
2 arcsin^/^ a 



ErT=l«nELo4-(CT) 



(28) 

(29) 

(30a) 
(30b) 
(30c) 
(30d) 
(30e) 

(31) 
(32) 
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Taking into account the form of the integrals in Eq. ((30)) we can express the deflection angle as 

(^i^ arcsin^/^ (cr) 



•\/ Fi (cr) + F2 arcsin a 



(33) 



2. Region 11 (a < z < 1) 

We now come to the second region where the potential is expressed in terms of a series around the local maximum 
at fj.. After expanding to first order, and using the same interpolating potential as in the case of the Schwarzchild 
metric Vb{z) = V{^) — p{z — /i)^, we obtain 



where 



and 



(1) 



dz 



1 



Vb{l)-Vi,{z) 



(34) 



(35) 



a = l — /X, b = z — n- 



(36) 



Since 



we have 



fe=0 



Afc(z) 



V{l)~V{z) 



Vb{l)-Vb{z) 
In terms of the new variable u = b/a one obtains 



(1) 



where = (/i — tr)/ (/_* — 1) and = 1. 
Notice that 



2 r+ du 



1 - 



-p(a + 6) 



Af,(z(w)) 



(a + &) 1 + M 



Therefore we get 



n=2 ' k=0 



where we have defined 



Jfc(T) 



- l(u+ 1) 



-du. 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 
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The first integrals are 



where 



The PMS yields 



and we obtain 



Jo(o- 
Ji{i 



e- 1 



e+ 1 
e- 1 



e + 1 
e - 1 



-ln(6) 



(e+ l)2e 



(l + 4e + e2)+ln(e) 



^^HH2,» + 14,^ + 2.-l)-|l„(.) 

(-'= + .» - 17.' - 62.' - 17.» + . - 1) + 1 1„(.), 



/i — (T + \/ a"^ — 2[LG + 2/i — 1 

^7^1 ■ 



-j^ cxD n— 1 

PPMS = 1^ " '^)""^ 



n=2 



21n3/2 



Once again, looking at the structure of the integrals (j43|) we write 



v/Gi(a)+G2 Ine 

Notice that the explicit expression of the coefhcients F\^2 and Gi^2 will depend on the metric. 



(43a) 
(43b) 
(43c) 
(43d) 
(43e) 

(44) 
(45) 
(46) 



(47) 



IV. RESULTS 



A. Schwarzschild metric 



Our first application is to the Schwarzschild metric, which corresponds to 



B{r) =A-'{r) = 1- 



(48) 



Here M is the Schwarzschild mass. The angle of deflection of a ray of light reaching a minimal distance rg from the 
black hole can be obtained using Eq. ([3]). The exact result can be expressed in terms of incomplete elliptic integrals 
of the first kind 31| and reads 



A</) = 4a 



where tq = tq/GM and 
T = 



/ro 


- 2 




+ 6 



, K = ^J{T -ro + 6)/2T , = / arc 



2 + T - r 



6 + T - r 



(49) 



(50) 



12 




FIG. 6: Left panel: deflection angle for the Schwarzschild metric as a function of /i using difi^erent approximations and the 
exact result. Right panel: percent error of the approximate solutions as a function of fi. (color online) 



[I] 



Method I 

Method I simplified 
Method II 



FIG. 7: Percent error for the deflection angle in the Schwarzschild metric calculated with Method I (Eq. (|25|l 'l. with the 
simplified version of Method I (Eq. ((53]) and with Method II (Eq. (|A8)) '). (color online) 



We compare our analytical formulas with the exact one and with the approximation proposed by Bozza 



2M- 



-2 In 



^(2+V3)(M 



i: 



(51) 



For brevity, we shall refer to the approximation developed in the previous Section as Method I. In the Appendix we 
have also derived a simpler analytical expression for the deflection angle using an alternative method, which we refer 
to as being Method II. 

In the former approach, the expression of Eq. (P5)l can be further simplified mantaining good accuracy by observing 
that 



2V6/i 



arcsm 2 



12 



/6arcsin(l-^);.2-8,.+ ^l^ 



-SW - 1 - 4 



(52a) 
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TABLE I: Numerical values of the coefficients A and B: ^Eiroa and BEiroa are taken from Table 1 of Ref. [lO|, whereas Aue 
and Bus are obtained from Eq. (|59|l . 

\q\ 0.1 M 0.25 M 0.5 M 0.75 M 1 M 

^Eiroa 2.00000 2.00224 2.01444 2.06586 2.19737 2.82843 

Ans 2.00000 2.00224 2.01444 2.06586 2.19737 2.82843 

BEiroa 0.207338 0.207979 0.21147 0.225997 0.262085 0.426782 

Bus 0.213892 0.214535 0.218032 0.232554 0.268419 0.430856 



and 

2^/6A^ liJ 



v/4Ai + l/31n -5— . (52b) 



Therefore we get the much simpler expression 

In Fig. [HI we compare our approximations, the one of Bozza |24|, and the exact result of Eq. (|15|) . It is clear that 
— apart from the region very close to fj. = 1 — our approximations are more accurate, especially as /i increases: 
the curves in the left plot corresponding to using our approximations are hardly distinguishable from the exact one. 
Indeed, in the case of Method I, after having reached a maximum error close to /i = 1, the precision keeps improving, 
as one can see in Fig. [71 In the case of Method II and of the simplified expression of Eq. ((53)) the error does not tend 
to zero as the photon sphere is approached but it appears to remain finite. 

B. Reissner- Nordstrom metric 

The Reissner-Nordstrom (RN) metric describes a black hole with charge and corresponds to 

i?(r)=A-i(r)= (^1-^ + ^) , D{r)^l. (54) 
The corresponding potential is found to be 



V{Z) = { trZ^ - 2 — Z + 1^ _ ^ «^^4^ (55) 



where rps — 4q^/(3A/ — ■\/9A/^ — 8q^) is the photon sphere and /i = ?'o/''ps; q is the charge of the black hole. Also 

in this case we have found that a ~ 1 — l/2/.t is a satisfactory analytical approximation to the optimal value upMS- 

Although the exact solution of the PMS condition would in general depend on the charge q, we have verified that 

the general features discussed in the case of the Schwarzchild metric remain valid and only a quite limited error is 

introduced by this choice. 

Using the former wc obtain 

-4^3 + - 4a + 2 -2a^ + - 6a^ + - Ga^ + Sa^ - 2a + 1 

Fx = 5 1 5 ^ (56a) 

3a(a2 + l)^ 12a(a2 + l)2 

ft^i + g, (56b) 
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FIG. 8: Percent error of the deflection angle in the Reissner-Nordstrom metric calculated with Method I at different values of 
/i. (color online) 



where a = ^/i/I" 



1, and 
Gi 

G2 



-3a^ + 6^2 + 1 (21a8 + \2of' - llOa^ - 28^2 - 23) g2 



3(a2-l)(a2 + l)' 
1 {Jo!^ - 10^2 + 31) 
4 



96 (a2 - 1) (a2 + 1)' 



128 



(57a) 
(57b) 



where the coefficients F\ 2 and G1.2 were previously introduced in the previous section (see Eqs. ((33|) and (|47|) ). 

Usin g ou r expression for the deflection angle we have obtained the coefficients A and B introduced by Eiroa ct al. 
in Ref. 



10|: 



where e = tq ,ps, , ps 
We have found 



A(/i K. -AXogiB e) - TT, 
rns, 7'ns bciug thc distance corresponding to the photon sphere. 



(58) 



A = - hm , ^f'PMS , 

l0g(/i- 1) 



cxp[ 



'2 /A/2 + 3v/9-8g2/M2 - 9 
0.234 + 0.203/V' - 1.096 vWW^+"(I826)], 



(59a) 
(59b) 



where V = ^9 - 8g2/Af2_ 

Table U shows that our analytical formulas are in remarkable agreement with the numerical results of Eiroa et 
al. [lol |. In particular our expression for the coefficient A appears to be exact. 

Fig. [8] displays the percent error of our approach for different values of /i as a function of q/M. Again, the error is 
generally below 0.5% and gets smaller as ji increases. 

Notice that we do not need to discuss the WDL of our formulas, since in this case one has that /i ^ 00 and the 
method of 26, 2^ is recovered. The reader will find a detailed comparison of our method with other methods available 



in thc literature in 



26 



27| 
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C. Janis-Newman-Winicour metric 



Finally, we consider the spherically symmetric metric solution to the Einstein massless scalar equations [32 1 
A{r) = {l-h/ry\ B{r) = {\-b/r)\ D{r) = (\ - b/rf'", 
which reduces to the Schwarzscliild metric for = 1 and h = GM . For this metric we obtain the potential 



V{z) = - 1 - 



bzV 



roj V ^0 



that can be expanded around z = to give 



(60) 



(61) 



V{z) « VlZ + V2Z'^ + v^z^ + . . . 

« _2(i._l)(l_l)2-il^ + 
ro ro 



1 - (:/- l){2iy-l) 1 



bV-Wb-^'' 



ro 



ro 



--(iy-l)iy(2iy-l){ — ] 1 

3^ ^ ^ ' \roJ \ ro ro 



z' + 0[z^] 



(62) 



Notice that the radius of convergence of the series of V{z) around z = is z = ro/b. We therefore ask that z > 1, i.e. 
that ro > b. Clearly the accuracy of the expansion above will depend upon the location of z and a larger number of 
terms is expected to be needed when z approaches one. 
We then find 

{a - 1) (a {?,Q^ + _ 23a + 16) - 13) 



(a-l)(l + a2) 



8a 



-V2 



64 (a^ + a) 



-V3 



F2 = '-l 



_ (6(a.-1)a.+ 1) ^, 

""1 — 7r~F7^ 7\ ^'3 



G2 



8^(2^ - 1) 

V2 + 3^W3 



(63a) 
(63b) 
(63c) 
(63d) 



where, again, a = — 1 and Fi^2 and Gi^2 have been introduced in Eqs. ([55]) and ([Tf|) . Notice that in the cubic 
approximation we can express the coefficient vi in terms of the other parameters: 



Vi = ~2flV2 -ifJ. V3. 



(64) 



Fig. [9] compares our approximation for the deflection angle with the exact result. One can see that also for this metric 
the accuracy is very good. 



V. CONCLUSIONS 



We have presented a new method for obtaining analytical expressions for the deflection angle of light in a static 
and spherically symmetric metric, which is accurate in both the weak and strong regimes. The former corresponding 
to lensing at arbitrarily large distances from the compact body and the latter to distances arbitrarily close to the 
photon sphere. Our first-order analytical formulas exhibit errors below 1% at all distances from the compact body. 
For this reason, although our method can be applied to any given order, in a way similar to what has been done 
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FIG. 9; Left panel: deflection angle in the Janis-Newman-Winicour metric using = 1/2 and & = 1; the solid line is the 
numerical result, whereas the dashed line corresponds to our analytical formula. Right panel: percent error for the deflection 
angle using our analytical formula. 



in Ref. [27| with the method of Ref. 
physical applications. 



26j, the accuracy of our first-order expressions is certainly sufficient for most 



26 



27| 



Moreover, the method that we have presented in this paper reduces to the previous method discussed in Refs 
in the weak deflection limit (WDL), since the arbitrary parameter a tends to unity in this limit. For this reason, for 
the comparison with alternative methods developed to describe the WDL we have relied on the discussion contained 



26 



27|. 



To the best of our knowledge, our method is the only one available that allows one to obtain completely analytical 
formulas which arc valid in both SDL and WDL regimes, regardless of the particular static and spherically symmetric 
metric used. The results that we obtain are clearly non~perturbative, since they do not correspond to a polyno- 
mial expression in any of the physical parameters in the model and provide the correct logarithmic strength of the 
singularities. Moreover, our analytical expressions never involve special functions and are easy to evaluate. Just to 
mention one success of our approach, in the case of the Reissner-Nordstrom metric, we have obtained an analytical 
formula for the coefficients A and B which have been numerically calculated by Eiroa et al. in Ref. 10[: in the case 



of the coefficient A our formula reproduces all the digits given by the numerical calculation. We also wish to mention 
that our method relies on solid mathematical grounds and could be used to obtain an exact series representation for 
the integrals of the deflection angle: we have not considered this issue because of the highly precise results that are 
obtained already working to order one. 

Finally, we have also discussed — showing an application to the Schwarzschild metric — an alternative method that 
yields larger errors but even simpler expressions and that can be useful in cases where a somewhat reduced accuracy 
can be traded for the possibility of more convenient analytical manipulations. 



APPENDIX A: ALTERNATIVE FORMULA THROUGH LINEAR INTERPOLATION 



In this appendix we present an alternative method to obtain the deflection angle in a given static and spherically 
symmetric metric. This method is a generalization of a recently published approach to the period of the simple 



17 

pendulum [33|. It is our purpose to obtain a simple analytical approximation to an integral of the form 

Jo VOW) 

where Q{1) = 0. We assume that any other zero of Q{z) is outside the closed interval [0, 1]. 
We define a reference function 

go(z) = (l-z)(z-zo), (A2) 
where zq < 0, and carry out the change of variables 



z = — — + — — sm 9. (A3) 



1=1 — 6^0 = arcsm 



We thus obtain 

We next substitute a linear funcion a + (36 for ^J F{z{9)) such that F(0) = a + /?6'o, F{1) = a + /?7r/2 and obtain 
the general approximate expression 



xjm\ 



(A5) 



We have yet to specify the exact location of zg- Notice that 



^a),i„ gill. (A6, 

(1 - z)(z- zo) zo - 1 

In the case of the Schwarzschild metric we have = y/QjII — tt, where / is given by Eq. (|Aip . with 

Q{z) = (Z - Zl)(z - Z2)(Z - Z3), 

zi < 0< Z3, Z2 = 1, (AT) 

and Z3 > 1 if /i > 1. 

We may set the location of zq to have either the most accurate analytical expression or the simplest one; in what 
follows we choose the latter. If zq = —1, then Oq = and 



« ^ ( + v/^^) In - n. (A8) 

Notice the logarithmic singularity at /i = 1 that comes from the fact that z^{ii = I) = Z2 = 1 and the integral diverges 
as /i ^ 1^. This approach is considerably less accurate than the preceding one, but we have decided to include it in 
this paper for two reasons: first, it provides simple and general expressions; second, its error is quite uniform. 
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